Flux tube in turbulent flow and quantum chromo dynamics 
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' Using some assumptions about the correlation function of velocity in a turbulent flow, a cylindri- 

cally symmetric tube-like solution is obtained. It is proposed that this turbulent flow is similar to a 
^jrj, flux tube in quantum chromodynamics. Using similar assumptions for quantum chromodynamics, 

r—{ • a flux tube filled with a color electric field is obtained. The connection between turbulence theory 

and quantum chromodynamics is discussed. 
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I. INTRODUCTION 



Two of the most serious problems in modern physics are the creation of turbulence and the quantization of strongly 
interacting quantum fields (nonperturbative quantization). From Ref. [l|: "In the world of physics, the difficulty of 
understanding turbulence in fluids is legendary. A famous quotation is also variously attributed to Einstein, Heisen- 
berg, Richard Fcynman. or Arnold Sommerfeld: Turbulence is the last great unsolved problem of classical physics." 
In a similar way the nonperturbative quantization is unsolved problems in quantum physics. The nonperturbative 
£Li quantization problem is connected with the problem of quark confinement or unobservability of a single quark. From 
Wikipedia 0: " Color confinement is the physics phenomenon that color charged particles (such as quarks) cannot be 
. isolated singularly, and therefore cannot be directly observed. The reasons for quark confinement are somewhat com- 
' plicated; there is no analytic proof that quantum chromodynamics should be confining, but intuitively, confinement 
is due to the force-carrying gluons having color charge. As any two electrically-charged particles separate, the electric 
O^l ■ fields between them diminish quickly, allowing (for example) electrons to become unbound from nuclei. However, as 
two quarks separate, the gluon fields form narrow tubes (or strings) of color charge, which tend to bring the quarks 
. together as though they were some kind of rubber band. This is quite different in behavior from electrical charge. 
C^ - ' Because of this behavior, the color force experienced by the quarks in the direction to hold them together, remains 
constant, regardless of their distance from each other." Here we would like to consider a possible connection between 
these two problems. 

In order to understand the reason of confinement problem in quantum chromodynamics there are many attempts to 
compare quantum chromodynamics witn other branches of physics. In Ref. Q it was suggested that the confinement of 
quarks into hadrons may happen due to a condensation of special gluonic configurations called Abelian monopoles. In 
this approach (the dual superconductivity scenario) a condensate of the monopoles breaks spontaneously an internal 
U(l) gauge symmetry. According to the dual superconductor idea, the breaking of the dual symmetry gives rise 
naturally to the dual Meissner effect, which insures a formation of a tube, which in turn leads to the confinement the 
quarks into hadronic bound states. In Ref. @ it is supposed that the SU(2) Yang-Mills theory in the low-temperature 
phase can be considered as a nematic liquid crystal. In Ref. [f| it is proposed that the SU(2) Yang-Mills theory can 
be interpreted as a two-band dual superconductor with an interband Josephson coupling. This approach is based on 
a slave-boson decomposition I5[ used in the high-T c superconductivity models (similar decomposition used in gauge 



theories is named as a spin-charge separation @ 0) 

The color "magnetic" flux tubes which appear in QCD with high baryon density in Ref's are considered. 

QCD with high density at low temperature is in the color-flavor locked(CFL) phase, where it exhibits color super- 
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conductivity. Unlike color electric/magnetic fluxes in gluon plasma, those color magnetic flux tubes are topologically 
and dynamically stable. They are a mixture of a superconductor vortex and a superfluid vortex and are called a 
semi-supcrfluid vortices. 

The color magnetic flux tubes also exist (as BPS vortices) in supersymmetric QCD. In Ref's [H[-[li[ one can nn d 
reviews on this subject. 



II. REYNOLDS NUMBER AND COUPLING CONSTANT 



Why such connection between turbulence and quantum chromodynamics may exists ? In both cases we have the 
fluctuations: velocity, pressure and so on for a turbulent fluid in the first case and fields (for example, SU(3) gauge 
fields for quantum chromodynamics) in the second case. In the turbulent fluid the fluctuations arc statistical ones 
and in quantum field theory the fluctuations are quantum ones. 

The character of fluid flow depends on Reynolds number Re 

Re = (1) 
M 

where p m is the fluid density, v is the fluid velocity, I is a characteristic length for a given flow and p is the fluid 
viscosity. If Re < Re cr then the flow is laminar one, if Re > Re cr then the flow is turbulent one. The open question 
in hydrodynamics is the creation of turbulence by Re rs Re cr . 

In quantum field theory there are perturbative regime when a dimensionless coupling constant 

o* = ±f (2) 
he 

is small enough a 2 < 1 and nonperturbative regime when a 2 > 1 (here g is a dimension coupling constant, H is Planck 
constant and c is the speed of light). In quantum field theory the open question is a nonperturbative quantization 
for a 2 > 1. For example, the fine-structure constant in quantum electrodynamics is a 2 = e 2 /hc w 1/137 (e is the 
electron charge), a 2 s» 0.1 for a weak interaction and a 2 > 1 for a strong interaction. 

Here we propose a possible connection between hydrodynamics and quantum field theory based on a comparison 
between Reynolds number (in hydrodynamics) and dimensionless coupling constant (in quantum field theory). Let 
us rewrite Eq. |T]) in following form 

Re = (3) 

in such a way that [/w 2 £ 4 ] = [l/.? 2 ] = g ■ cm 3 /s 2 , [h] = [pi 3 ] = g ■ cm 3 /s. Now we will assume that there is following 
relation 

1/g 2 <- pv 2 l\ (4) 

h <-> pi 3 . (5) 

It means that we have following correspondence 

• h = <-» p = 0. In this case a classical theory corresponds to an ideal fluid. 

• h ^ and a 2 < 1 <-> p ^ and Re < Re cr . In this case a laminar fluid corresponds to a perturbative regime of 
quantum field theory 

• h ^ and a 2 > 1 <-> p ^ and Re > Rc cr . In this case a turbulent fluid corresponds to a nonperturbative 
regime of quantum field theory 



III. FLUX TUBE IN TURBULENT FLOW 



How we can use proposed here the similarity between hydrodynamics ans nonperturbative quantum field theory ? 
One can try to find in hydrodynamics an analog of flux tube from quantum chromodynamics. 

At first we will try to examine a static ideal and laminar fluids. The Navier - Stokes equations are 

dv 

— + {vV)v = Vp + vAv (6) 

at p m 
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here we consider incompressible fluid divv = and v — is the kinematic viscosity. If the viscosity p = then we 
have the ideal fluid. We want to find cylindrically symmetric solution of the Navier - Stokes equations but not in a 
pipe. For the ideal [v = 0) and laminar fluids the equation is 



dv, d v 



Vp—r- = v 



dp dp 1 



(7) 



but v p — for the flow along axes z. Consequently the solution for the ideal fluid does not exist. For the laminar fluid 
we have d 2 v z /dp 2 = and this equation has a solution in a pipe only. Only one possibility exists to have a nontrivial 
solution: it is necessary to have the correlation between v p and dv z /dp in the LHS of Eq. ([7|. 

Now we would like to consider a turbulent fluid. We will try to find an axially symmetric solutions. We consider 
the fluid in unbounded space with two sources located at the z = ±oo. Then the equations (J6j> are 

Id, . dv z dv z dv z . . 

v - p T P {pVz) = v '-te +v >-at + v *W {) 

Id dv p dvp dv p 

%T P {pVp) = v '-te +v >-at + v *W {) 

We assume that in the first approximation such flow can be described with fluctuating velocities v zp that depend on 
p only. We use the cylindrical coordinate system z, r, <fi. The flow is along the axes z: v z ^ and v p = 0, (■ ■ ■ ) means 
the statistical overaging. Remember that we consider the solutions in unbounded space, not in a pipe. 
Now we can write equation for v z velocity 



Id, dv 

'-pT P [pVz) = Vp l-p 

Unfortunately it is not too easy to write equation for v p velocity. For example, if we write equation in the form 

Id, „ dv 



%t P ^ )= <w (12) 



and after averaging wc will have the correlation 



dv 

v p ^=0 (13) 

because v p = 0. It happens because the terms v z ^- and v^^- are not taking into account in Eq. ©. 

In order to write equation for v p we multiply Eq. ([9]) on v p and only after that we average over the statistical 
fluctuations. In this case wc assume that one can neglect with the terms V z ^- and w</>^- 

vv p v" p + -v p v' p = v 2 p v' p (14) 

where (■■■)' — d(- ■ ■ )/dp. After some manipulations and averaging wc have 

2 VP j''-^W+^W= 1 rH (15) 

here we neglect with d/dz and d/d(j> terms in Eq. ©. Now we have to do some assumptions about correlations 
functions in (|15p . Wc should have the information about 2 and 3 points correlation functions. The most important 
simplification is that all correlation ructions are described by one function ^(x 1 ), i = 1,2, 3: 



- V 



v a {x*)v b {t) = S ab ^(x l )^(y l ), (16) 



dv^z 1 ) ( ■ dcj)(z l ) 
v a (x l )v b (y % ) c = lA a bc4>(x l )(f>(y 1 )—^- h transmutation by (x,y, z) 

(a a c/>(j/ I )(/)(z l ) + transmutation by (x,y, z)) + 

f3 a (f)(y l ) — — h transmutation by (x, y, z) ) (17) 

oz l J 
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here A abc and a a , Pa are some constants, a,b,c = p, <f>. The correlation function (| 1 T[) is choosen in such way because we 
assume that v a = and conjectural corrections to the first term 4>(x 1, )4>(y l ) ^ can be only quadratic compositions: 
either or </>(?/') d< ^ z i ^ . For v z we have to use not v z but: v z = v z — v z 



vaix^Mv 1 ) = i<t>{x l )4>{f), (18) 



Va(x i )v b {y i )v z (z l ) = 7 [{A abc (j)(x l )(l){y l )(l){z' 1 ) + transmutation by (x, y, z)) + 

(aQ(/)(?/ l )0(z l ) + transmutation by (x, y, z))] (19) 

where 7 is a constant indicating the difference in the scale between longitudinal velocity v z and the fluctuations of 
transversal velocity v p . 
According to (JTHJ) dTTJ) 



v 2 p = 2 , (20) 
3</>V + 2a<j) 2 + 2/3#' '. (22) 



Thus we have following equation 



<lV = #' + ^ + k (23) 
p 3 3 

Now we would like to write equation for v z = v z + Sv z velocity component. After substitution into Eq. ((lip , averaging 
over statistic fluctuation wc have 



. --v P Sv' z (24) 

here the RHS is similar to the RHS of Eq. (fT2")l and consequently similar calculations leads to 

vv" z + - p v' z = 1 Utf + \a<$> + l^'J (25) 

here the number 7 appears in the consequence of correlations (|T5|) p^)) . Thus we have following set of equations for 
the longitudinal velocity component v z and fluctuating transversal velocity component v p 

4>" + = W + a^ + j3(t)', (26) 

v" + ^v' = y fa' + &4> + P4>') (27) 

where we have introduced dimensionless variables x = p/po, po = v/4>07 <i> = 'fi/'t'o, S = ^av/4>o, (3 = |/3/<^ , 
v = v z /(j)Q and ^>o = 0(0). The numerical solution of the set ([26]) (|27|) is presented in Fig. [1] Let us note that the 
solution exists only for negative values a < 0, j3 < 0. 

Our interpretation of this solution is follows. There exists a turbulent flow moving in z direction. In the flow exists 
a core (tube) with the velocity profile given in Fig. [TJ curve 1; the z— component of the velocity at the infinity is 
i'oo ^ 0,p — > 00; the profile of statistical fluctuations describing by the function <f> is given by curve 2 in Fig. [2 as 
4> — > then at the infinity the flow become laminar one. ft means that we have to choose the solution in such a way 
that Woo becomes small enough that flow becomes laminar one. In our opinion it is an analog of a flux tube filled with 
a gauge fields and stretched between quark and antiquark from quantum chromodynamics. 

Such kind of solution is characterized by having a strongly pronounced core (tube) through which flows the turbulent 
fluid. Such flux of the turbulent fluid may be inserted either into a flow of a laminar fluid Vrx, ^ (that is more 
probably) or in a rest fluid = (that is less probably). 

In summary we list the basic ideas and assumptions in presented approximation: 

• All degrees of freedom are splitted into two parts: ordered movement v z and disordered fluctuations v p . 

• The correlation functions of the second and third orders are expressed in terms of one function (one-function 
approximation) . 

• The correlation function of the third order is expressed in terms of the correlation function of the second order. 
The last assumption allow us to cut off an infinite equations set relating all correlation functions of all orders. 
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FIG. 1: The profiles of functions v(x) and 4>(x), a = j3 = —3,7 = 6. v(x) - curve 1, (p(x) - curve 2 

IV. FLUX TUBE IN QUANTUM CHROMODYNAMICS 

In this section we would like to show that similar calculations can be done in quantum chromodynamics. The basic 
assumptions will be same as listed at the end of Section IIIII 

In this section we follow the conventions of Rcf. Q . Starting with the SU (TV) gauge group with generators T B we 
define the SU(N) gauge fields, A^ = A B T B . Let G be a subgroup of SU(N) and SU(N)/G is a cosct. Then the 
gauge field A^ can be decomposed as 

= A B T B = a°T a + A™T m , (28) 
a° G G and A™ G SU(N)/G (29) 

where the indices a,b,c... belongs to the subgroup G and m, n, . . .to the coset SU(N)/G; B are SU(N) indices. 
Based on this the field strength can be decomposed as 

rf v T B = T« u T a + T™T m (30) 

where 

>%« = d »< - d >< + f abc ^< G G, (32) 

= f amn A™A n v G G, (33) 

Kv = F ?u + G™ G SU(N)/G, (34) 

F™ = d^A™ - d u A™ + r np A^Al G SU{N)/G, (35) 

G™ = f mnb (A»a b v - Ayl) g ST/(iV)/G (36) 

where f ABC are the structural constants of SU(N). The SU(N) Yang-Mills Lagrangian can be decomposed as 

£■ -F hv-Fbiiv — J~i lv J~im V + J~ UV ^F m(j,v (37) 

We will specialize to: 

• the SU(2) case when we let SU(N) -> 5*7(2), G -> C/(l), and / ABC -> e ABC ; the indices a = 3, m,n = 1,2 
(here e^ 5 * 7 = 0, ±1 is the absolutely antisymmetric Levi-Civita tensor). 
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• the SU(3) case when we let SU(N) -> SU(3), G -> SU(2); the indices a = 1, 2, 3,m, n = 4, 5, 6, 7, 8. 

For the first case one can assume that the fourth order Green function can be decomposed as the product of the 
second order Green functions and obtain Ginzburg-Landau Lagrangian [l6| . 

For our purposes the second case is more interesting: using all assumptions from Section HTTl with some modifications 
we will show that in SU(3) gauge theory exists a flux tube. In this case a longitudinal color electric field E\ will be 
analog of the velocity v z and gauge potentials A™ in the coset SU(3)/SU(2) are the analog of fluctuating transversal 
velocity v p . 

In quantizing the classical field equations system (via Heisenberg's non-perturbative method [I3|) one first replaces 
the classical fields by field operators A B — ► A B . This yields the following Yang-Mills equations for the operators 

d u P Bflv = 0. (38) 

These nonlinear equations for the field operators of the nonlinear quantum fields can be used to determine expectation 
values for the field operators A B , where (• ■ • ) = (Q\ ■ ■ ■ \Q) and \Q) is some quantum state. One can also use these 
equations to determine the expectation values of operators that are built up from the fundamental operators A B . For 

example, the "electric" field operator, Ef = d Af — d z A B + gf BCD A^Ai 7 giving the expectation (£ B ). The simple 
gauge field expectation values, (A^x)), are obtained by average Eq. (|38[) over some quantum state \Q) 



Q 



d„£ BfM/ 



Q) = 0. (39) 



One problem in using these equations to obtain expectation values like (A B ), is that these equations involve not only 
powers or derivatives of (A B ) (i.e. terms like d a (A B ) or d a dp(A B )) but also contain terms like Q B ^ = (A B A^). 
Starting with Eq. (f3"9"| one can generate an operator differential equation for the product A B A^ thus allowing the 
determination of the Green function Q B ^ 

(Q \A B {x)d yv T B ^{x)\ q) = 0. (40) 

However this equation will in turn contain other, higher order Green functions. Repeating these steps leads to an 
infinite set of equations connecting Green functions of ever increasing order. This construction, leading to an infinite 
set of coupled, differential equations, does not have an exact, analytical solution and so must be handled using some 
approximation. The problem is the same as in the turbulent theory. 

For obtaining a tube-like solution in SU(3) gauge theory we do similar assumptions as in Section [TTT] (here we follow 
to |l|): 

• The gauge field components A® G SU{2) belonging to the small subgroup SU(2) are in an ordered phase: 

(a;(x)) = a%{x). (41) 

oJJ is the analog of v z . 

• The gauge field components A™ £ SU(3)/SU(2)) are in a disordered phase: 

(A™(x))=0. (42) 

• The 2-point Green function 6*2(2;, y) approximately are 

G 2 (x, y) = (A™(x)A?(y)) = -^/ rapo / n **^ o (x)0 i (») > p,v = 1, 2, 3, (43) 
here <p a is a real triplet scalar fields. They are the analog of fluctuating velocity v p . 

• Schcmaticaly the 4-points Green function G&{x, y, z, u) can be decomposed by following way 

G 4 (x,y,z,u) = ^A™{x)Al{y)Al{z)A q f3 {u)^ = \G 2 {x,y)G 2 {z,u) + permutations by (x, y, z, u)] + 

[S mn r] l j /V G2(z, u) + permutations by (x, y, z, u)] . (44) 
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After some calculations one can reduce initial SU (3) Lagrangian to 

C eff = -- (T^T A » V ) = - \f; v F^ + \ (D^ a ) (D^ a ) - ^AlA c ^ b - 

where a, A, m, (M 2 )^ v are some parameters. The details can be found in Ref. [HI]. Corresponding field equations are 

D u F afiv = ge abc 4> b D^4> c - (M 2 )^ A a u , (46) 
D^4> a = -\<j) a {<t> b 4> b - /Lt 2 ) - uA\ A a ^cj) b (47) 

where = + ge^A 1 ^ is the covariant derivative. Using following ansatz for Eq's (|4"6")l (|4"7| 



A(p) = ^; Aj{p) = ^; <f(p) 

9 9 9 (48) 

(M 2 Y V =diag{M o ,M 1 ,0,0} 



gives us 









h — 




z 








z 








z 



'/' + ^ = « (0 2 - f - M 2 ) , (49) 
/ (0 2 + 1; 2 - M 2 ) , (50) 

6" + ^ = 0[-/ 2 + « 2 + A(0 2 -m 2 )]. (51) 

The numerical investigation of the field equations shows that there exists a regular cylindrically symmetric solution 
with nonzero color electric fields 

El = ~fv, El =--/'• (52) 
9 9 

The flux of the chromoelectric field E\ is nonzero 

oo oo 

$ £ = 2tt J E 3 z pdp = 2tt y /upd/o 7^ (53) 



which is similar to nonzero flux of a turbulent fluid across a core in Section IIIII 



27r / p m (v z - Woo) pdp ^ 0. (54) 



o 

V. CONCLUSIONS AND DISCUSSION 

From author point of view the problems of turbulence and quantum chromodynamics are connected with an infinite 
equations set connecting all correlation (for turbulence) or Green (for quantum chromodynamics) functions. This 
equations set can not be solved exactly. Heisenberg wrote that there is only one way to solve such equations 
system: to cut off the set using some assumption on a Green function and to decompose some Green function as the 
combination of Green functions of lower orders. This is exactly what we have presented in this paper. 

Now we would like to enumerate common features being used here for finding tube-like solutions presented above: 

• For cutting off an infinite equations set for correlation/Green functions the correlation/Green function of 
third/fourth order is decomposed via lower correlation/ Green functions. 

• All correlation functions are expressed via one function. The Green functions are expressed via two set of scalar 
fields. 

It is necessary to emphasize that between the turbulence and quantum chromodynamics there is a difference. Any field 
theory has Lagrangian and field equations are corresponding Eulcr equations for this Lagrangian. The Navier-Stokcs 
equtions are not Euler equations in the consequence of the presence of viscocity. 
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